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Abstract 

A general discussion is made concerning the ways in which one can get signatures 
about a possible liquid-gas phase transition in nuclear matter. Microcanonical tem¬ 
perature, heat capacity and second order derivative of the entropy versus energy 
formulas have been deduced in a general case. These formulas are exact, simply ap¬ 
plicable and do not depend on any model assumption. Therefore, they are suitable 
to be applied on experimental data. The formulas are tested in various situations. 
It is evidenced that when the freeze-out constraint is of fluctuating volume type the 
deduced (heat capacity and second order derivative of the entropy versus energy) 
formulas will prompt the spinodal region through specific signals. Finally, the same 
microcanonical formulas are deduced for the case when an incomplete number of 
fragments per event are available. These formulas could overcome the freeze-out 
backtracking deficiencies. 
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1 Introduction 


From a long time it is believed that, due to the van der Waals nature of the nucleon-nucleon 
interaction, nuclear matter is likely to exhibit a liquid-gas phase transition [1], While such a 
phase transition was predicted within various models, from the experimental point of view 
so far things are rather inconclusive. The reasons are related to both the inherent difficulties 
of backtracking the freeze-out stage and the lack of any rigorous mathematical apparatus 
permitting the unambiguous determination of a liquid-gas phase transition when it exists 
in experimental data. Some years ago, the general conviction was that the experimental 
measurement of the nuclear caloric curve could give a def ini te answer to this problem. Ne¬ 
glecting the weaknesses of each of the “thermometers” tested with that occasion, there is one 
fundamental barrier preventing the inference of a first order phase transition from any exper¬ 
imentally obtained caloric curve even with back-bending shape: The experimentally obtained 
microcanonical systems follow an unknown path in the excitation energy ( E )- freeze-out vol¬ 
ume ( V ) plane [2], Even if the experimental path intersects the coexistence region the phase 
transition may remain unrevealed in the caloric curve [2]. For overcoming that situation it 




was proposed that the heat capacity being a quantity independent on the experimental path 
but depending on the local freeze-out constrain could give precise information concerning 
the occurrence of the transition [3,2], More exactly, if the experimental path is intersecting 
the spinodal region of the system’s phase diagram, this would be signaled by a negative 
value of the heat capacity. An analytical heat capacity formula was proposed with that oc¬ 
casion and was subsequently employed in investigating the experimental data obtained by 
the MULTICS-MINIBALL collaboration occasion with which negative branches of the heat 
capacity curves were identified [4]. The formula was further applied on the INDRA data with 
the same success [5]. However, it is worth noticing that though formally correctly deduced, 
this formula is based on unevaluated terms (such as the microcanonical temperature (T) or 
the heat capacity (Ci) corresponding to a partial energy of the system (E\)) which further 
imply a chain of assumptions and approximations Q resulting in (unpredicted) deviations 
from the exact results. 

The above matters motivated us to draw the attention in the present paper on an alternate 
way to calculate the microcanonical heat capacity of nuclear systems proposed in Ref. [6]. As 
will be further shown this method is exact, model independent, simply applicable and does 
not depend on any unevaluated term. The present paper is structured as follows: Section 2 
gives a brief review of the model used for illustrating the paper’s ideas. Section 3 investi¬ 
gates on the possible signals one can (experimentally) get about a possible liquid-gas phase 
transition. Microcanonical formulas for temperature, heat capacity and second-order deriva¬ 
tive of the entropy versus energy are deduced for various conservation options in section 4. 
Examples for the functioning of the above-mentioned formulas are given in section 5. The 
same quantities are deduced in section 6 for the case in which only a percent of the total 
number of fragments are available per event. Conclusions are drawn in section 7. 


2 Model review 


For testing the paper’s argumentation we use the sharp microcanonical multifragmentation 
model proposed in Ref. [7]. The model concerns the disassembly of a statistically equilibrated 
nuclear source (A, Z, E , V) (i.e. the source is defined by the parameters: mass number, atomic 
number, excitation energy, and freeze-out volume respectively); its basic hypothesis is equal 
probability between all configurations C : {A*, Z % , e*, r*, p*, i = 1 (the mass 

number, the atomic number, the excitation energy, the position and the momentum of each 
fragment i of the configuration C, composed of N fragments) which are subject to stan¬ 
dard microcanonical constraints: J2i A* = A, )T, ; Z, = Z, Pi = 0, r « x p* = 0, E - 
constant: integration over the fragments’ momenta can be analytically performed in the ex¬ 
pression of the total number of states, then one works in the smaller configuration space: 
C' : {Aj, Z % , €i, rj, i — 1,..., iV}; a Metropolis-type simulation is employed for determining 
the average value of any system observable. Function of the desired conservation restriction 
(TV) level (n) one can write: 


1 For example in Ref. [4] the temperature estimator does not correspond to the microcanonical 
one, is based on the assumption of thermal equilibrium between fragment internal and external 
degrees of freedom and, moreover, depends on additional level density parameters; furthermore, C\ 
is deduced from the T((E\)) dependence which results in even larger deviations. 
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with n = 0,1,2; ft(0) = S(H-E), 11(1) = 11(0) <$(EiPi), ft(2) = ft(l) <5(E^ x Pi ); 
/(0) = 1,/(1) = l/(Ei '^i) 3 ^ 2 and /(2) = /(1)/Vdetl, where I is the inertial tensor of the 
system. Here H denotes the system’s Hamiltonian: H = E?; P 2 /(2m*) +Ei<j hq+E?: Q — E?: Pi 
and K = E — J2i<j Vq — E?; Q + Ei -B? (Tq is the Coulomb interaction between fragments i and 
j; and are respectively the internal excitation and the binding energy of fragment i). 
In the standard version of the model fragments are idealized as nonoverlapping hard spheres 
placed into a spherical recipient; intersection between fragments and the “recipient’s” wall 
is also forbidden. We call this freeze-out hypothesis (i). In a simplified version, (ii), one 
can further perform an approximate integration over the position variables as in [8] and get 
the extra factor: Vf ree = fl;=i (V — i V 0 /N) in the expression of the statistical weight of a 
configuration C" : {A i: Zj, e,, i = 1, (i.e. Wc» = WcVf ree ), where V 0 stands for 

the volume of the nuclear system at normal nuclear matter density. In this version of the 
model the Wigner Seitz approximation is used for the interfragment Coulomb interaction. 
When angular momentum conservation is also considered, det I is approximated by the 
quantity corresponding of a sphere of volume V uniformly filled with nuclear matter of 
density (V 0 /V)p 0 . The latter version is computationally faster and provides a different freeze- 
out perspective (i.e. spherical hard-core interaction is missing). 


3 Echos of the liquid-gas phase transition 


Which are the echos one can experimentally get from a liquid-gas phase transition? For an¬ 
swering this question one should first define the liquid-gas phase transition in finite extensive 
or nonextensive systems. This matter was recently addressed in [8]. There, it was shown that 
the techniques usually employed for very large systems work for the case of finite extensive 
or nonextensive systems as well. Namely, it was shown that curves like Y(X)\y* (here X 
is an extensive variable, Y is the X’s conjugate and Y ! is another intensive variable) have 
the property of revealing the phase transition (even in small and nonextensive systems) by 
bending backwards. The microcanonical isobaric caloric curves or isothermal pressure versus 
volume curves are natural examples for the physical situation under study. The probability 
distributions of an isobaric canonical ensemble proved to be a precious tool for identifying 
the phase transition when it exists (by exhibiting a double peaked structure) and for con¬ 
structing the corresponding T(E)\p and P(V)\t curves. Reasons lay in the simple expression 
of the probability of a state with excitation energy E and volume V in a constant pressure 
canonical ensemble characterized by the parameters P (canonical pressure) and /3 (inverse 
of the canonical temperature): 
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leading to the expression of the microcanonical entropy (S = In W(E, V)): 
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where Z can (/3, P ) represents the constant pressure canonical partition function. Starting from 
their definitions (T /4 = (dS(E,V)/dE)~ 1 , P ^ = T li dS(E,V)/dV), one can easily obtain the 
microcanonical temperature and pressure: 
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One can then follow the method proposed in [8] and follow a P^ = P or a = 1/P 
path in the ( E , V) plane by using the (P ;j , Tfl) pairs evaluated (by means of eq. (4)) in each 
(E, V ) point with sufficiently collected events from the probability distribution of a canonical 
ensemble at constant pressure (one can simulate this ensemble within the present model by 
multiplying the statistical weight of a configuration by e~ /3E ~ l3PV and letting both energy 
and freeze-out volume to fluctuate). The corresponding isobaric and isothermal paths are 
plotted in Fig. 1 (upper-right panel) for a system with A = 200, Z = 82 for which Coulomb 
interaction is switched to zero. The corresponding constant pressure canonical probability 
distribution is represented by a contour plot (in log scale) in the same picture. The double 
peaked structure points the presence of a first order phase transition. Both isothermal and 
isobaric paths are intersecting the two peaks of the distribution and also a third saddle-like 
point. The corresponding T(E)\p and P(V)\t curves are plotted in the lateral panels with 
thick lines. They both present a backbending region indicating the phase transition. 


While such constructions can be simply performed in theoretical models, experimentally 
one can usually only access microcanonical ensembles (due to the sorting of events in bins 
of E) . Moreover, it is widely believed that the events corresponding to a microcanonical 
ensemble at a given E have a fluctuating volume at freeze-out - somewhat similar to the 
case of a volume constraint e -AV proposed in Ref. [2] (where A is the Lagrange multiplier 
corresponding to the volume observable). As clearly resulting from Fig. 1, the constant A 
microcanonical ensemble corresponds to the folding on the V axis of the constant pressure 
canonical ensemble. Indeed, replacing in eq. (2) pP with A, and making the integration over 
V one gets: 
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and then the expression of the microcanonical temperature: 
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The dependence of the last two quantities on E is represented in Fig. 1 (lower panel). Thus, 
V X (E) is a projection on the E axis of the V(E, V ) distribution. Therefore, a double peaked 
structure in V X (E) is necessarily related to a double peaked structure in V(E , V ), and thus, 
is reflecting a first order liquid-gas phase transition as defined in Refs. [9,8]. 


Of course, the same reasoning may be done for the projection on the V axis. The correspond¬ 
ing curves are plotted in Fig. 1 (left panel). However, the relevance of this latter projection 
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Fig. 1. Right-upper panel: probability distribution (in log scale) of a canonical ensemble at con¬ 
stant pressure corresponding to the nucleus (200,82) where Coulomb interaction has been switched 
off. The canonical parameters /3 and P are labeled on the picture. The thick lines correspond to 
microcanonical isobaric and isothermal paths corresponding to P^ = P and T ^ = 1//3. The lateral 
panels correspond to projections (dashed lines) of the probability distributions on V/Vq (left) and 
E (down) axes. Canonical pressure (left panel) and constant-A temperature are figured with thin 
lines. The corresponding microcanonical curves are plotted with thick lines. (Here and in the rest 
of the figures E stands for excitation energy.) 

for the experimental data is conditioned by a good event by event resolution in freeze-out 
volume which was not achieved so far. 


4 Microcanonical formulas 


Let us further deduce the exact expressions for microcanonical temperature and heat capac¬ 
ity. The microcanonical density of states corresponding to a total energy of the system E 
writes: 
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Here F(V) is the factor in front of f{n) from eq. (1) (V is a generic notation for fragment 
partition) and £ = J2i<j Vij + the expression (9) being formal. Then, eq. (7) 

may be rewritten as: 
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The micro canonical temperature writes T — dS/dE with S = In W(E). Which means: 
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If the limits of the (formal) integral over £ are not depending on E, then one can simply 
make the derivative versus E inside the integral from eq. (10): 


T~ i 


1 

W(E) 


£w c , 

C' 


/ |(JV-n)-l \ 


^ |(JV-n)-l \ 


( 12 ) 


where we used the implication: K = E — £ and the notation( ) for the average over the 
ensemble’s states. The heat capacity of the system is by definition: C 1 = — T 2 (d 2 S/dE 2 ). 
Using the same arguments as those from the deduction of T one obtains: 
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Alternatively, one can (similarly) evaluate the second order derivative of the system’s entropy 
versus E: 
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In these terms, the system’s spinodal region will be pointed by positive values of the above 
quantity. One can easily check that for n — 1 (i.e. energy and total momentum are the 
conserved quantities), eq. (14) is similar with the one deduced in Ref. [6]. The microcanonical 
formulas given by eqs. (12), (14) and (15) are universally applicable to any system for which 
the “external” degrees of freedom can be treated classically irrespective to other specificities 
of the system. Indeed, the integral over the momenta variables (with the various possible 
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Fig. 2. Comparison between the results of the microcanonical formulas (for temperature, heat 
capacity and second order derivative of the entropy) and the corresponding curves deduced from 
the probability distributions of a canonical ensemble with a given freeze-out constraint (i.e. constant 
volume or constant A). Each column correspond to the case labeled on it. 

conservation options) given in eq. (1), from which the microcanonical T, C and d 2 S/dE 2 
formulas were further deduced is valid for any classical N particle system. This is the case 
of nuclear multifragmentation as well: it is widely accepted that, since at freeze-out the 
system is rather dilute, a classical treatment of the fragments’ motion is appropriate. These 
features make eqs. (12), (14) and (15) suitable to be applied for experimental data: they 
are independent on specific model assumptions (i.e. freeze-out hypothesis, internal excitation 
treatment - level density, interfragment interactions, etc.), and moreover, they only depend 
on two parameters (K and N) which have to be estimated for the freeze-out stage. 


5 Examples 


While the deduction of the above formulas guarantees their accuracy, some examples of their 
functioning on concrete cases are further discussed. As an independent test, the results will 
be compared with the temperature/heat capacity computed using the method synthesized by 
eq. (6). More exactly, the caloric test curves are computed using eq. (6) from the probability 
distributions of a canonical ensemble with a given freeze-out constraint and the heat capacity 
and d 2 S/dE 2 curves are farther evaluated from the obtained caloric curves (C = dE/dT, 
d 2 S/dE 2 = (—l/T 2 )dT/dE). Calculations are performed with both versions of the model 
((i) and (ii)), for two different freeze-out constraints: constant A and constant volume, for 
the case in which the Coulomb interaction is present and the one in which it is switched 





K(M eV/u) 


Fig. 3. Kinetic energy distribution corresponding to the source nucleus (200,82) with the freeze-out 
volume V = 10 Vo with the Coulomb interaction included, for three excitation energies (0.5, 1 and 
2 MeV/u). 

off. (This variety of examples constitute a further test for the universality of the deduced 
formulas.) The considered conservation “level” is n = 2, i.e. the total energy, the momentum 
and the angular momentum are conserved quantities. Results corresponding to the cases: 
(50,23) nucleus with constant A = 5 10 -4 fm -3 with Coulomb interaction included, version 
(ii) of the model, (200,82), A = 3 10” 3 fm -3 , without Coulomb interaction, version (ii) of 
the model and the same nucleus with Coulomb interaction included, but a constant volume 
constraint: V — 10 Vo, version (i) of the model are presented in Fig. 2 (first, middle and 
third column respectively). In the first two cases, backbendings of the caloric curves reflected 
in negative branches of the heat capacity curves and positive regions in the second order 
derivative of the entropy versus E - all related with a first order phase transition can be 
observed. In the third case (constant volume constraint) the caloric curve has a monotonic 
increase with a “plateau-like” region (where the curve slope is smaller, but positive) which 
is reflected in a positive C (with a peak in the plateau-like region) and a negative d 2 S/dE 2 . 
(As will be further seen, one cannot draw any conclusion from the third case concerning 
the intersection of the coexistence region since a constant volume path is not necessarily the 
order parameter of the system.) A very good agreement between the curves calculated using 
the probability distribution of a canonical ensemble and the points evaluated by means of 
the microcanonical formulas is to be noticed. The calculated microcanonical heat capacity 
and second derivative of the entropy versus excitation energy points are clearly revealing 
the spinodal region. While the calculated heat capacity points have large error bars near the 
border of the spinodal region due to the asymptotic behavior of the heat capacity, d 2 S/dE 2 
has small error bars even in the border region - thus being an even more robust quantity for 
identifying the spinodal region. 

Some details concerning the applicability of the microcanonical formulas have to mentioned. 
As mentioned earlier, the validity of these formulas is restricted by the condition that the 
integration domains in eq. (10) over the quantity £ not to be restricted by the total energy 
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Fig. 4. Caloric curves, heat capacity curves and second order derivative of the entropy curves 
evaluated with the microcanonical formulas for the case of the nucleus (200,82), without Coulomb 
interaction and a freeze-out constrain corresponding to a path in A given by: A (E) = A ^(E, 10 Vo). 


of the system E. One should thus limit the cases on which this formula is applied to those 
obeying the above condition. But how to distinguish those “forbidden'' points from the rest? 
A restriction by E of the upper limit of the integral over £ has to be clearly pointed by a “cut” 
in the right part of £ distribution. Since K = E — £, a similar “cut” (corresponding to K = 0) 
has to be noticed in the left part of the K probability distribution. In other words, when 
the probability that K = 0 is significant (i.e. a cut clearly affects the K distribution) then 
deviations from the microcanonical formulas should occur. These “cuts” (usually) occur at 
very low energies, so in a region unimportant for the system’s thermodynamics. All the points 
represented in Fig. 2 have been checked against the cut in the K probability distribution. For 
illustrating this idea we take the case (200,82), with Coulomb, V — 10 Vo. The probability 
distributions of K corresponding to the E =0.5, 1 and 2 MeV/u cases are illustrated in 
Fig. 3. While the K probability distributions corresponding to E =0.5 and 1 MeV/u are 
intersecting the K = 0 axis, (giving thus a non-negligible probability for K = 0), the 
probability distribution corresponding to the E = 2 MeV/u case is not intersecting the 
K = 0 vertical axis. This fixes the lowest excitation energy limit for which the microcanonical 
formula works to some value between 1 and 2 MeV/u. 

The crossing of the spinodal region by a (experimental) path in the (E, V) plane may not be 
reflected as a backbending in a caloric curve. However, if the volume constraint is of the e~ AV 
type the spinodal region is likely to be revealed by a negative branch in the heat capacity. 
This point was nicely argued in Ref. [2] and can be evidenced by means of the present 






microcanonical formulas as a further test for their abilities. To this aim, we chose A as the 
microcanonical “A” parameter (A M ) corresponding to a constant volume path {V = 10 Vo) 
in the (E,V) plane. I.e.: 

X(E) = A jL1 (E, 10 Vq). (16) 

The microcanonical A can be simply obtained from eq. (4): 
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As illustrated in Fig. 4 this path appears to intersect the system’s spinodal region without 
inducing any backbending in the caloric curve for the case of the system (200,82) without 
Coulomb interaction, (ii) version of the model. Indeed, while the caloric curve has a mono- 
tonical increase, the corresponding heat capacity presents a negative branch and the second 
order derivative of the entropy a positive region. The negative region in the heat capacity and 
the positive region in the second order derivative of the entropy versus energy are evidencing 
a “spinodal line” of the system corresponding to the volume V — 10 Vq in the system’s 
phase diagram in (E, V ) representation. The limits E of this region are approximately those 
delimited by the above mentioned signals. 



6 Fewer fragments 


Let’s now assume that experimentally only a (small) number N x (< N) of fragments cor¬ 
responding to a given event are detected. Can one still get the desired thermodynamical 
information about the system? This question is addressed in the present section. Let us turn 
back to eq. (1). For the highest conservation level (n = 2) we can write: 
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Pi — Ei^l r i X p,;. The integration over the momenta of the first Ni fragments can be further 
performed: 
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where Ii is the inertial tensor corresponding to the first Ni fragments at freeze-out, M x = 
rrii and 0 stands for the step function. One can further deduce the thermodynamical 



quantities obtained before similarly, starting from their statistical definitions. In the general 
case (i.e a generic “level” of conservation n) one obtains: 
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where Jl 2 = K x - JC(n) with /C(0) = 0, /C( 1) = Pj/{2M 1 ) and /C(2) = /C(l) + Li T Ii 1 L 1 /2. 
Following the same argumentation as in section 4 one can deduce the criterion of validity 
of the above formulas: the probability distribution of K 2 must not intersect the K 2 = 0 
vertical axis at a non-negligible value. The advantages of the above formulas are obvious: 
they only depend on one parameter (K 2 ) and, more importantly, they can be used for 
inferring information about the liquid-gas phase transition even when one has incomplete 
information about the fragmentation events (i.e. only N { fragments from the N fragments 
of one fragmentation event are detected). For a given microcanonical formula, depending on 
the desired conservation level (n ), the minimum allowed value of N\ is hxed by the condition 
that the factors containing N\ be positive. Thus, for example, when only the total energy, 
E, is conserved (n = 0) one only needs one freeze-out fragment per event in order to obtain 
the correct microcanonical temperature! 


7 Summary 


In summary, the possibilities of obtaining information about a possible liquid-gas phase 
transition from experimental heavy ion collision data have been investigated. The connec¬ 
tion between negative heat capacity (or positive dS 2 /dE 2 ) obtained in a fluctuating volume 
ensemble and the liquid gas phase transition has been discussed. It was shown that a negative 
value of the heat capacity in a fluctuating volume (constant A) ensemble is in connection 
with a first order phase transition as defined in Ref. [8]. Microcanonical formulas for tem¬ 
perature, heat capacity and second order derivative of the entropy versus energy have been 
rigorously deduced for various conservation options (conserved energy; conserved energy and 
momentum; conserved energy, momentum and angular momentum). These formulas were 
tested with an independent method for obtaining the above mentioned quantities which 
extracts the microcanonical information from the probability distributions of a canonical en¬ 
semble with a given freeze-out volume constraint. An excellent agreement between the two 
methods have been obtained for all the considered curves ( T(E ), C(E), d 2 S/dE 2 (E)) for 
various cases ((50,23), A = 5 10~ 4 frn -3 with Coulomb, (200,82), A = 3 10~ 3 frn -3 , without 
Coulomb, (200,82), V — 10 Vo with Coulomb). In the first two cases, backbendings in the 
caloric curves reflected in negative branches of the heat capacity curves and positive regions 



in the d 1 2 3 4 5 6 7 8 9 S/dE 2 (E) curves have been obtained reflecting the occurrence of a first order phase 
transition. In the third case no backbending in the caloric curve or negative region in the 
heat capacity curve was observed which however does not imply any conclusion about the 
occurrence of a first order phase transition simply because a constant volume path in the 
(E, V ) plane is not necessarily the system’s order parameter. Indeed, choosing for the nu¬ 
cleus (200,82) without Coulomb a dependence of A on E corresponding to the microcanonical 
lambda with V — 10 Vo (A (E) = A ^(E, 10 Vo)) hi an ensemble with fluctuating volume one 
gets a well defined negative branch in the heat capacity curve and a positive region of the 
second order derivative of the entropy in spite of a monotonically increasing caloric curve. 
The case of incomplete (experimental) information concerning the fragmentation event is 
further discussed. Microcanonical formulas corresponding to all conservation levels depend¬ 
ing only on a given number of fragments (smaller then the total number of fragments from 
the given event) are deduced. While demanding some more kinematical information about 
the detected fragments (only for the n > 0 cases) their utility is obvious. Thus, for example 
in the n — 0 case one only needs the kinetic energy of one freeze-out fragment per event in 
order to deduce the correct microcanonical temperature. The resulted microcanonical for¬ 
mulas are model independent thus being applicable on experimental data. Moreover they do 
not depend on any unevaluated term (such as heat capacity of a partial energy or the level 
density parameters as is the case in Ref. [4]). Of course, the success of these formulas depends 
on the accuracy of backtracking the information about the primary fragments. Here the mi¬ 
crocanonical formulas based on partial break-up information may be very effective since 
some particles (such as freeze-out neutrons) carry at the asymptotic stage the untouched 
information about the freeze-out. 
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